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The Ekpyrotic scenario is studied in the context of some extensions of Hofava-Lifshitz gravity. 
Some particular solutions that lead to cyclic Hubble parameters are analyzed, where the corre- 
sponding gravitational actions are reconstructed by using several techniques and auxiliary fields. Its 
comparison with standard F(R) gravity is performed. In addition, the so-called Little Rip, a stage 
of the universe evolution when some bounded systems may be dissolute, is also studied in this frame 
of theories. 
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I. INTRODUCTION 

During the last years, one of the main theoretical problems that has concerned to the scientific community refers to 
the mystery of the accelerating expansion of the universe. Since a deviation in the luminosity distance of Supernovae 
la was observed in 1999, along with other independent observations (such as Cosmic Microwave Background (CMB) 
anisotropics), the fact that the expansion is accelerating has been mostly accepted. Such acceleration seems to be 
approximately an effect of a cosmological constant in the Einstein field equations, exactly the same form that the 
vacuum energy density acquires. Nevertheless, the observed value is so small in comparison with the one predicted 
by quantum field theories that the probable existence of a dark component in the universe has been established, 
which would be responsible for the accelerating expansion, and may deviate in principle from a perfect de Sitter 
acceleration (in other words from a cosmological constant). Under the name of dark energy, plenty of candidates 
have been proposed which can perfectly predict the observational data, leading to a problem of degeneracy. On 
the other hand, standard model of cosmology also requires an initial stage, the so-called inflation, in order to solve 
some cosmological problems as the homogeneity or the flatness problems. As an alternative, the so-called ekpyrotic 
scenario may avoid the need to provide initial conditions (inherent in every inflationary model), since the universe 
evolution acquires a periodic behavior, such that in every cycle a new universe is born (see Ref. [l|). In addition, it 
is argued that the problem of flatness does not appear in this model because the universe initially was in a nearly 
BPS (Bogolmonyi-Prasad-Sommerfield) state, which is homogeneous (see Ref. Q). In the last years, very promising 
models capable to unify the entire cosmic evolution under the same mechanism have been proposed, where the 
inflationary epoch and the late-time acceleration era are unified under the same mechanism (or alternatively the 
ekpyrotic scenario), providing a simpler picture of the universe evolution. Most such models are described by scalar 
fields due to its simple form (see Ref. 0] and references therein), or other kind of fields (see Ref. 0), but also a large 
effort has been done in the reconstruction of modified gravity theories (for a general report, see Ref. 0]) available to 
reproduce the cosmic evolution (for a review, see Ref. [|, and Refs. HQ), which 

may seem more natural as they are 

expressed in terms purely of the metric tensor without additional fields. 

On the other hand, a new theory of gravity that is power-counting renormalizable has been proposed recently in 
Ref. Q. Such theory, already known as Hofava-Lifshitz gravity, breaks the invariance under full diffeomorphisms of 
General Relativity by introducing an anisotropy between the spatial and time coordinates through a critical exponent 
z. This restriction of the symmetries allows the theory to be power-counting renormalizable, but an additional scalar 
degree of freedom is found, which introduces instabilities in the spectrum of the theory (see Refs. [l,[l(|). However, 
some extensions of the theory seem to address the problem of the scalar mode [TJ E3 > as we ^ as *° generalize the 
action to more complex ones (see Ref. [IH). Moreover, cosmological models have been widely studied in the context 
of Hofava-Lifshitz gravity (see Ref. and also generalizations of the original action (similarly to standard F(R) 

gravity) have been proposed, where the entire cosmological history can be well reproduced, and it has also a good 
UV behavior (see Refs. [H, El). 

The aim of the present paper is to study the ekpyrotic scenario in the frame of some extensions of Hofava-Lifshitz 
gravity, where a universe described in terms purely of gravity is able to pass along the different stages of an ekpyrotic 
model. This class of cosmological solutions can be realized in standard F(R) gravity as shown in Ref. 17:]. Here, we 
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reconstruct some periodic solutions for the Hubble parameter, which may be able to describe the entire evolution of 
the universe. In addition, we also analyze the shape of the action for each phase of the ekpyrotic scenario, where 
the possibility of the occurrence of a Little Rip is explored. The so-called Little Rip is a postulated phase of the 
universe evolution, when a very strong accelerating expansion would lead to break some bounded systems, as the 
Solar System or even the molecules and atoms (see Ref. [Uj]). Such breaking is shown to be fully compatible with 
the ekpyrotic scenario in comparison with future singularities as the Big Rip that are not, unless some cure for the 
future singularity is considered [19(. Moreover, the presence of a Big Bang/Crunch singularity, usual in ekpyrotic 
cosmologies, is still an open issue for this kind of cyclic scenario, where quantum effects may resolve it (see [20|) or 
an effective theory that generates a non singular bounce (see [2l|). Nevertheless, here we are interested to explore 
the classical effects of the theory, where some non singular solutions are proposed, while the study of possible UV 
effects in the presence of the singularity is beyond the purpose of this paper. 

The paper is organized as follows: in the next section, F(R) Hofava-Lifshitz gravity is briefly reviewed. In section 
III, the actions for some cyclic solutions are reconstructed. Finally, section IV is devoted to the analysis of ekpyrotic 
scenario, where each phase of the cycle is analyzed. 



II. MODIFIED F(R) HORAVA-LIFSHITZ GRAVITY 

In this section, modified Hofava-Lifshitz F(R) gravity is briefly reviewed [l3|, [TH, [l6[ . We start by writing a general 
metric in the so-called Arnowitt-Deser-Misner (ADM) decomposition in a 3 + 1 spacetime (for more details see |22j|). 

ds 2 = -N 2 dt 2 + g$ (dx l + N l dt)(dx 3 + N 3 dt) , (1) 

where i,j = 1, 2, 3, N is the so-called lapse variable, and N l is the shift 3-vector. In standard general relativity (GR), 
the Ricci scalar can be written in terms of this metric, and yields 

R = K l3 K lj - K 2 + R^ + 2V Al KV„n 1/ - n»V„n>*) , (2) 

here K = Kij, Kij is the extrinsic curvature, R^ 3 ' is the spatial scalar curvature, and a unit vector perpendicular 
to a hypersurface of constant time. The extrinsic curvature Kij is defined as 

Ki^T^idlf-V^Nj-VfNi) . (3) 

In the original model Q , the lapse variable N is taken to be just time-dependent, so that the projectability condition 
holds and by using the foliation- preserving diffeomorphisms ([5]), it can be fixed to be N = 1. As pointed out in 
imposing the projectability condition may cause problems with Newton's law in the Hofava gravity. For the non- 
projectable case, the Newton law could be restored (while keeping stability) by the "healthy" extension of the original 
Hofava gravity of Ref. [ll[ . 

The action for standard F(R) gravity can be written as 



S = J d A xy/g&)NF{R) . (4) 

Gravity of Ref. Q is assumed to have different scaling properties of the space and time coordinates 

x l = bx\ t = b*t , (5) 

where z is a dynamical critical exponent that renders the theory renormalizable for z = 3 in 3 + 1 spacetime dimensions 
0. GR is recovered when z = 1. The scaling properties (0 render the theory invariant only under the so-called 
foliation-preserving diffeomorphisms: 

Sx* = C(x\t), St = f(t). (6) 

It has been pointed that, in the IR limit, the additional scalar degree of freedom can be removed by means of an 
additional U(l) symmetry pjj]. Here, we are interested on actions as follow, 

5* = i / dtd 3 xy^)NF(R) , R = KijK ij - XK 2 + R^ + 2/xV M (fi"V ' v n v - n v V v n^) - L (3) (ffg } ) , (7) 
2k J 
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where k is the dimensionless gravitational coupling, and where, two new constants A and /x appear, which account 
for the violation of the full diffeomorphism transformations. Note that in the original Hofava gravity theory [8|, the 
fourth term in the expression for R can be omitted, as it becomes a total derivative. This generalization of the Hofava- 
Lifshitz action, similar to standard F(R) gravity, may provide the way to describe the entire cosmological evolution 
with no need to introduce any additional field but where an additional scalar mode is assumed. The possibility of 
violations of Newtonian law, due to the extra scalar mode coming from F(R), can be avoided by the appropriate 
expression for the action, as it was pointed out in Ref. [l6|. In addition, standard F(R) gravity (U]) can be recovered 
by setting A = fi = 1. The term L^(g^) in the action ([7J is chosen to be Q 

L( 3 \gl?)=E«G ijk iE kl , (8) 
where the generalized De Witt metric is given by, 

In Ref. [1], the expression for Eij is constructed to satisfy the "detailed balance principle" in order to restrict the 
number of free parameters of the theory, and it is defined through the variation of an action 

Vj®Efi = , (10) 

°9ij 

The action Wfgjy] is assumed to be defined by the metric and the covariant derivatives on the three-dimensional 
hypersurface J2t- I n 0> ^[si?] i s explicitly given for the case z = 2, 

W = ±- fd 3 x^(R-2A w ), (11) 
K w J 

and for the case z = 3, 



W = \ f lo 3 (T) . (12) 



Here kw in (1111) is a coupling constant of dimension —1/2 and w 2 in (|12l) is the dimensionless coupling constant. 
W3(r) in (IT21 is given by 

„, ( r, = Tr (r a dr + a r a r) = £ * (rffvt. + ?r 3 r' m r £ ) a». . («) 

Here we are interested in the study of cosmological solutions for the theory described by action ([7]). Spatially-flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) metric is assumed 

3 

ds 2 = -N 2 dt 2 + a 2 (t)J2(d xl ) 2 , ( 14 ) 
»=i 

where N is taken to be just time-dependent (projectability condition) and, by using the foliation-preserving dif- 
feomorphisms (JB]), it can be set to unity, N — 1. Then, just as an assumption of the solution, N is taken to be 
unity. 

For a flat FLRW metric ([H]) , and a vanishing cosmological constant, the scalar R is given by 

3(1 - 3A + 6/i)ff 2 6^ d (H 



+ ^UJ- (15) 

For the action ([7]), and assuming the FLRW metric (TT51) . the second FLRW equation can be obtained by varying the 

(3) 

action with respect to the spatial metric , what yields 

= F{R) - 2(1 - 3A + 3/i) [h + 3i? 2 ) F'{R) ~ 2(1 - 3X)HRF"{R) + 2fi (r 2 F (3 \R) + RF"{R)^j + n 2 p m , (16) 
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here k 2 = 16irG, p m is the pressure of a perfect fluid that fills the universe, and N = 1 . Note that this equation turns 
out the usual second FLRW equation for standard F(R) gravity (QJ when A = fj, = 1. If we assume the projectability 
condition, variation over N of the action ([7j yields the following global constraint 







d 6 x 



F(R) - 6F'{R) |(1 - 3A + 3p)H 2 + fillX + 6fiHRF"(R) 



(17) 



Now, by using the ordinary conservation equation for the matter fluid p r , 
it yields 



'3H(p m +p m ) = 0, and integrating Eq. ([16 



= F(R)-6 (1 -3\ + 3p)H 2 + pH F'(R) + 6pHRF"(R) 



K 2 pr, 



«3 



(18) 



where C is an integration constant, taken to be zero, according to the constraint equation (|17|) . In [23], however, 
it has been claimed that C needs not always vanish in a local region, since (JXTj) needs to be satisfied in the whole 
universe. In the region C > 0, the Ca~ 3 term in (|T5|) may be regarded as dark matter. 

If we do not assume the projectability condition, we can directly obtain (fT8|) . which corresponds to the first FLRW 
equation, by varying the action ([7]) over N. Hence, starting from a given F(R) function, and solving Eqs. (TIB)) and 
([17}. a cosmological solution can be obtained. 



III. RECONSTRUCTING CYCLIC UNIVERSES 

The aim of this section is to show that any cosmology may be realized in F(R) Hofava-Lifshitz gravity. For this 
purpose, we present two different methods of reconstruction, the first one is based on the use of the number of 
e-foldings, while, the second one uses an auxiliary scalar field. 



A. Reconstructing cyclic universe using e-folding 



We will assume the flat FLRW metric defined in (|14p with N = 1, in such a case the first FLRW equation is given 
by (fT8|) with C = 0, which can be rewritten as a function of the number of e-foldings rj = In ^ instead of the time t. 
This technique has been developed in [24j for classical F(R) gravity, and for Hofava-Lifshitz F(R)-gravity Since 



% = H4- and #7 = H 2 -^ 

at clt] ar) z ar\ A 



H^-f, the first FLRW equation (22) is rewritten as 



= F(R) - 6 



-H 2 + fiHH' 
3 



dF{ ^ ' 6fiH 2 [2AHH 1 + 6nH' 2 + ^H"H'} - p , 



dR 



d 2 R 



(19) 



where A = 3 — 9A + 18/i and the primes denote derivatives respect 77. By using the energy conservation equation 
p + 3H(1 + w)p — 0, the energy density yields, 



p = poa- 3 ^ = p ao 3il+w) e- 3 ^+ w ^ 



(20) 



As the Hubble parameter can be written as a function of the number of e-foldings, H — H(rj), the scalar curvature in 
(fl~5l) takes the form 



R = AH 2 + QpHH' 



(21) 



which can be solved respect to 77 as rj — rj(R). Then, the equation ([T^]) for F(R) with the variable R is obtained. 
This can be a little simplified by writing G(rj) — H 2 instead of the Hubble parameter. In such a case, the differential 
equation (JT9J) gives 



= F(R) - 6 



dR 



6/1 [AGG' + 3/jGG" 



d 2 R 



p a- 3(1+w) e^ 1+w ^ 



(22) 



and the scalar curvature is now written as R = AG + 3/iG". Hence, for a given cosmological solution H 2 = G(rj), one 
can solve the equation ([2"2"1) . and the corresponding F(R) is obtained. 
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In order to illustrate that cyclic solutions can be reproduced by this kind of theories, let us consider the following 
example: 



ff(*) = -^xsinfe 



(23) 



where Fx and T are constants. The number of e-foldings is: 



„, . Ida 2tt „ . /2tt \ da 2tt tt . /2tt , , 
H(t) = -- = -—Fx sm -t =► - = --Fx sm ] A 



77(i) = In 



a(t) 
a 



Hi 



2tt , 
cos [ —t ]-l 



Using (jlll A[) . the function G(r?) and its derivatives are given by: 



(24) 



G(ry) = if 2 
Then, we have: 



(2Fx + ry) 7? , G'fa) 



,//,-„). (,'",,„ 



(25) 



ft = -3 ( ^ ) [2/iffi + 277 (/i + (1 - 3A + 6^) Fx) + (1 - 3A + 6/x) r? 2 ] 



Fx ± 



1-3A + 6/1 V \ (l-3A + 6/i) 



f 2 - 



R 



3(1 - 3A + 6/i) (^) 



(26) 



Now, if we call x = ±. jj 



Ht 



(X-3A+6^ ^ "1 3(1-3A+6 M )(^)" 



1 - 3A + 6/i 



we can write: 



Fx + x. 



(27) 



We also have that: 



dF(R) _ 1 dFi(x) 

dR 6(1-3A + 6m)( 2 t) 2 ^ 



d 2 F(R) 
dR 2 



6(l-3A + 6/x) (^) 2 a; 



d 2 Fx(a;) l<Fi(x) 
dx 2 x dx 



(28) 



where i?i(x) = F(R(x)). 

We can now rewrite (|25j) in terms of the new variable x by using (|27|) . leading to: 



Gfo(z)) = - I ^ 



T) V (1 - 3A + 6//)- 



~ Hi — 



2 V , 2 

-X + x 



1 - 3A + 6/x 



G'( V (x)) = 2 [% 



T ) V 1 - 3A + 6a* 



(29) 
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Finally, by introducing QUI AH29]) into the equation (1221) and considering the case of vacuum, we arrive to the following 
differential equation for F\ (x) : 



= x 2 F 1 (x) + 



1 — 3A + 6ju 



Hl- 



(f -3A + 6/i)< 



ly 2 



(f - 3A + 6p) 2 



H 1 )x - x i 



dFijx) 
dx 



1 - 3A + 6/i 



2 M , 2 

X + X 



(f - 3A + 6/Li) 2 J 1 - 3A + 6fi 



d 2 F 1 (x) 
dx 2 



(30) 



Here, we have obtained an equation for the gravitational action, that in principle can not provide an exact expression, 
but which can be integrated numerically. Hence, this solution reproduces a periodic behavior for the Hubble parameter 
leading to a cyclic universe. 



B. Reconstructing cyclic universe using a scalar field 



In this subsection it will be shown how to construct an F(R) Hofava-Lifshitz gravity model realizing any given 
cosmology, this time using instead the technique of @. We start from the action for F(R) Hofava-Lifshitz gravity 



S 



which is equivalent to 



5' 



dtd 3 xV^)N{F(R) + Cmatter), 



dtd?x^g&N(P{4,)R + Q{<j>) + C matter ) 



(31) 



(32) 



Here, C ma tter is the matter Lagrangian density and P and Q are proper functions of the scalar field, <fi, which can 
be regarded as an auxiliary field, because there is no kinetic term depending on <f> in the Lagrangian. By varying the 
action with respect to tf>, it follows that 



Q = P'(<f>)R + Q'(<f>), 

which can be solved in terms of <p, as 

^ = </>(&)■ 

By substituting ([34]) into (|32|) and comparing with (f3~T|) . one obtains 



(33) 
(34) 



S = / dtd 3 xy/^N{F(R) + C matter ) ? 



F(R) = P(cf)(R))R + Q(<f>(R)). 



(35) 



We proceed now in the same way that we did in Section II, assuming the FLRWL metric, the second FLRWL equation 
can be obtained by varying the action (|32l) with respect to the spatial metric g\j . This equation can be written as: 



P(4>) |E - 2 (1 - 3A + 3^) (ZH 2 + H^j | - 2 (1 - 3A) H 



dP{4>) d 2 P{(t)) 



dt 



2fi- 



dt 2 



Q(0)+P = O 



(36) 



If we assume now the projectability condition, we can obtain a global constraint doing the variation of the action ([7]) 
over TV, it yields: 



P{<t>) 1 7? - 6 (1 - 3A + 3/i) H 2 +fiH | + 6vH^j^ + Q{<j>) -p = 
We can combine (I36|) and (|37[) in order to eliminate the function Q(4>), we finally obtain: 

2 M ^» _ 2 (1 - 3A + 3 M ) H d -^1 -2(1- 3A) HP W t)) +P + P = 



dt 



(37) 



(38) 



7 



As we may redefine the scalar field <p properly, we can choose 

(j) = t. (39) 

Provided the scale factor a is given by a proper function g(t) as 

a = a e s{t \ (40) 

with a constant ao, and if it is moreover assumed that p and p are the sum of the different matter contributions, 
with constant equation of state (EoS) parameters w;, Eq. (|38|) then reduces to the following second order differential 
equation 

2M^^ - 2 (1 - 3A + 3a*) - 2 (1 - 3A) . 9 "(0)P(0) + £(1 + ^)p l0 a^ 1+ ^ e~^+^ = (41) 

From this equation we can obtain P(4>) and using Eq. (I37[) we find that 

Q(0 = -P(^) {ii - 6 [(1 - 3A + 3 M ) P 2 + M iPj } - + ^ paa^ 1 ^- 3 ^'^ (42) 

i 

As a result, any given cosmology, expressed as (|40|) . can indeed be realized (as anticipated) by some specific /(P)- 
gravity. Note that Eq. (|4*Tj) is a second order differential equation on P((j>) when g'((j)) is known, but it can also be 
considered as a first order differential equation on g'((j>) (i.e. on H(<fi)) in the case that the function P(<fi) is given. In 
the following we will use this last point of view to find out a function P(P) that reproduces a cyclic universe. 
When matter can be neglected Eq. (|4"Tj) can be rewritten as: 

±( mm ^)=^ m ^*m ,43, 

which can be solved as Q: 



/i 1 dP(0) 3^ 2 1-3A+3 M 



, Pi,,) ^ / ri*P(<^ > (^^Mj (44) 



1 - 3A P{4>) d4> (1 - 3A) 
In the second equality, we have used the partial integration. Furthermore by writing P(0) as: 



P((f>) = U(4>) i-^+tl (45) 



()III B| is rewritten as follows: 



, U s = 2p 1 dU(<P) 12/i 2 1 f,,(mS\ f4fil 

3W l-3A + 3/xP(0) # (1 - 3A + 3^) 2 U{<j)) 2 J 9 \ deb J ' 1 ' 

We now consider the case given by: 

2(1-3X1 2(1-3A) 

P((j)) = U{(f>) i-3A+3 f . = p [cos(oj(t>)} (47) 
where Pq and w are constants. Then, using Eq. (|43[) . the solution is given by: 

g'(4>) = 90 [cosM)] 2 + i _ ^ tan(w^) (l - 1 _ ^ + 3 ^ [sin^)] 2 ) (48) 

where go is an integration constant. Note that the tangent term in (|48|) makes the solutions to contain some divergences 
that correspond to points where the scale factor becomes null, i.e. a(to) = 0. These divergences can be identified 
with a Big Bang/Crunch singularity and they are very common in cyclic universes, where the ekpyrotic scenario is 
reproduced. In order to have a smooth transition through the Big Bang/Crunch singularity, one expects that the 
quantum effects of the theory will avoid the occurrence of the singularity. However, this is a large task, even more in 
a background solution as (14"8"1) . and should be explored separately in the future. In addition, other mechanisms for a 
smooth transition have been suggested as the introduction of an additional term in the action or a different coupling 
with the matter lagrangian (see Ref. (25|). 
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IV. EKPYROTIC SCENARIO IN HORAVA-LIFSHITZ GRAVITY 



We have shown above that periodic solutions can be easily reconstructed in the frame of extended Hof ava-Lifshitz 
gravity. Here we are more interested to analyze ekpyrotic models in such kind of theories. The so-called Ekpy- 
rotic/cyclic universe is an alternative explanation to the inflationary paradigm proposed one decade ago in Ref. [f|, 
that can provide a realistic picture of the universe evolution (for a confrontation between both models, see [Hj]). In 
the same way as the inflationary scenario, ekpyrotic cosmological models can also predict the origin of primordial 
inhomogeineties that leads to the formation of large structures and the anisotropics observed in the CMB. In addition, 
this model does not require initial conditions in comparison with the standard inflationary scenario due to its cyclic 
nature. In general, the cosmological evolution presented by an ekpyrotic universe consist of infinite cycles, where each 
cycle contains four stages: a first initial hot state similar to the standard Big Bang model, then a phase of accelerated 
expansion, after which the universe starts to contract and finally the cycle ends in a Big Bang/Crunch transition, 
when the cycle starts again. The cosmological problems enumerated above are solved during the contracting phase. In 
the usual ekpyrotic models, brane scenarios or scalar fields are considered (see However, it is clear that modified 
gravity, and precisely F(R) gravity, can perfectly reproduce the ekpyrotic scenario [l7|. Here we are interested to see 
how the cosmological problems can be solved during the contracting phase in the context of Horf ava-Lifshitz gravity, 
and to reconstruct the corresponding behavior of the action during each phase of an ekpyrotic universe. The first 
FLRW equation is given by, 

3 tt2 _ 1 ( PmO PrO AxO _ k \ 

where the subscripts refers to matter (m), radiation (r), anisotropics (a), and k is the spatial curvature, while Pf(r) 
is defined as, 

= w^^m) & {k) ~ ' i * HF ' lff> + 3 " H kF " l * > ) ■ (50) 

In order to solve the initial cosmological problems, the last term in ([4"!Jf should dominate over the rest when the scale 
factor tends to zero, i.e. when the universe approaches the Big Bang (Crunch) singularity. Hence, the effective energy 
density defined in (1501) should behave as Ppim °c l/a m with m > 6 when the scale factor tends to zero, such that 
close to the initial singularity, the FLRW equation (|49l) can be approximated as, 

3 C 

^H 2 ~ PF(k) ~ > (51) 

where C is a constant. Then, we can reconstruct the form of the action F(R) close to the Big Bang (Crunch) 
singularity by solving the FLRW equation. Hence, for the Hubble parameter (jSTj) . the scalar curvature is given by, 

R = [(I - 3A + 6/V) - ^m]n 2 — . (52) 

And the FLRW equation (|5T|) yields an expression where F(R) is the unknown quantity, 

k2F , l{k) + 2^(1 -3A + 3/1)-^^ _ {l-3\ + 6ri-ym R = Q 

2/1771 2/J777 

This is an Euler equation that can be easily solved, and gives the function for F(R), 

F(R) = k x RP+ + k 2 R p ~ . (54) 

where. 



3m/j - 2k 2 {1 - 3A + 3/x) ± J4n 2 (l - 3A + 3/j)(k 2 (I - 3A + 3u) - 3mu) + mp(8 - 24A + (48 + m)n) 
P± = 4^ • (55) 

Note that the scalar curvature tends to infinity when a — > 0, and in such strong gravity regime, the parameters A and 
fi should be different than one, the limit of General Relativity, as the breaking of Lorentz invariance will be present 
in such kind of regimes, while it is recovered for the weak field systems. Moreover, in order to get a smooth transition 
along the singularity, the first derivative of F(R) should tend to infinity to ensure that the matter energy densities 
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remain finite in (|49|) . which can be easily achieved when ((3± — 1) < in (|53|) . 

After this contracting phase, the ekpyrotic model suggests that a hot initial state, similar to the Big Bang model, is 
created (in the original ekpyrotic model by the collision between branes), and which may be created by the decaying 
of the extra scalar modes coming from F(R) in this class of theories. Nevertheless, this is beyond the purpose of this 
paper, where our aim is to show the approximated form that the action should look like for each phase of the cycle. 
Then, during the matter/radiation dominated epochs, the action may seem as the standard Hilbert-Einstein action 
with F(R) ~ R and R — R, i. e. the parameters responsible of the breaking of full diffeomorphisms should recover 
the values of GR, A = fJ, ~ 1. The last phase for each cycle refers to an accelerating era, which may be described by 
the usual ACDM model, whose Hubble parameter can be written in terms of the number of e-foldings as, 



H 2 = Hi 



-p a a 



-p % 3 e 3 " 



(56) 



where Ho and po are constants. In the frame of General Relativity, the terms in the r.h.s of equation (|56|) correspond to 
an effective cosmological constant A = 3iJg and to a pressureless fluid. The corresponding F(R) can be reconstructed 
by following the steps described above. For this case the function G(rj) is given by 



Giri) = Hi + y/>oa - 3 e" 3 " . (57) 
And by using the expression for the scalar curvature R = AG + 3/iG", the relation between R and rj is obtained, 



-37, 



R - AHg 



(58) 



where k — 
expression, 

1 + 3(/i - A) 
6/x(l - 3A) 



—poa ■ 



fc(3 + 9(/x-A)) ' 

Then, by substituting (|57|) and ([58]) in the equation ([22|) . one gets the following differential 



F(R) 



l + 3(/z-A) ~ 
3/x(l - 3A) 



3H 2 p{l - 3A + 6/x) 



2/x(l - 3A) 



F'(R)- {R-9pHl){R-3Hl{l-3\ + 6p))F"{R) = , 

(59) 



here we have neglected the contribution of matter for simplicity. By performing a change of variable x — 



d 2 F dF 
= i(l - x)-— + (7 - {a + /3 + 1) x) — 
dx A dx 



R~9fJ,Hg 



the equation (1591) can be easily identified as an hypergeometric differential equation, 

a/3F , 

Ul" UX' 

with the set of parameters (a, /3, 7) given by 



3ffg(l+3(M-A))' 



■P = 



1 + A(9/i - I) 



a/3 = - 



1 + 3(m - A) 
6^(1 - 3A) 



' 2(1 + 3(/x - A)) ' " ' ^ 3p(l -3A) 
The solution of the equation (|60p is a Gauss' hypergeometric function [l6j . 

Fix) = CiFia, (3, 7; x) + C 2 x 1 " r F{a - 7 + 1, j3 - 7 + 1, 2 - 7; x) 



(60) 



(61) 



(62) 



where G\ and C2 are constants. Then, this action reproduces the ACDM model described by the Hubble parameter 
(|56|) without including a cosmological constants. Note that this is the same result obtained in [24| for classical F(R) 
gravity, although in this case the solution depends on the parameters of the theory (/1, A) whose values differ from the 
classical theory. 

Other kind of accelerating expansions can be also reconstructed in the context of this class of theories as showed in 
Ref. However, due to the periodic behavior of ekpyrotic universes, models containing future singularities (usually 
phantom models) are not allowed in this kind of models unless a mechanism for avoiding the singularity is introduced. 
Nevertheless, a new class of phantom models that do not contain Big Rip singularities but only affects to bound 
systems without reaching a singular point, the so-called Little Rip, has been proposed in Ref. and extended to 
modified gravities in Ref. 17J . Basically, these cosmological models consist on a phantom-like evolution, free of future 
singularities but whose strong expansion breaks the bond of some coupling systems (as galaxies, solar systems, or even 
atoms, nuclei...), what has been called as a Little Rip. An simple example of this kind of evolution can be described 
by the Hubble parameter, 



Hit) ~ H t , 



(63) 
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where Hq is a constant. In this case, we can also reconstruct the corresponding F(R) action by solving the FLRW 
equation (fTS)) . The scalar curvature is given by, 

R = 3(1 - 3\ + 6fx)H^t 2 +Q/j,H . (64) 

Then, the FLRW flj§ yields, 



^FiR) - ( 3H ^ + (1 3A + y^ ® H ov) \ F , {R) + 6HQfi{R _ q Ho ^ f » {r) = o . ( 65) 



2 v ' V 1 - 3A + 6// 

This is also an hypergeometric equation, whose solution is given by, 

F(R) = [c 1 C/( 7 ,/3;x( J R))+C 2 4 Q )( a ;(^))] (R-6H Qf i)^ 2 , (66) 

where U('j,/3;x) is the confluent hypergeometric function and lSf\x) is the Laguerre polynomial. The variable x(R) 
and the set of parameters (7, /3, a) are defined as, 

_ (l-3\ + 3n)(R-6H ») _ (2 - 6A + 3 M ) fl _5 _3 
Xl j 6JJ M(l-3A + 6/i) ' 7 2(1-3A + 3 M )' P 2 ' * 2' [ ' 

Hence, the F(R) action (|66|) corresponds to a series of powers in R that are capable to reproduce a kind of behavior 
given by the Hubble parameter (j6"3")l . In such case, we have that the effective energy density can be approximated as, 

p F(R) (xt 2 . (68) 

Note that for a cyclic universe, as the ones studied in section above, the phase when the universe expansion is 
accelerated can be approximated by (|63[) . such that a Little Rip may occur in the ekpyrotic scenario. In order to 
show in a qualitative way how this Little Rip occurs, i.e. how some bounded systems are broken, let us compare the 
effective energy density ()68|) with the energy density of some known systems as the Solar-Earth system, and calculate 
the time remaining before the Little Rip occurs. By assuming that AF(fl)fto) = -tt-Hq ~ 10~ 47 GeV 4 , where the age 
of the universe is taken to be to ~ 13.73 Gyrs, according to Ref. [27j, and a mean density of the Sun-Earth system 
given by p©-© = 0.594 x 10~ 3 kg/m 3 ~ 10~ 21 GeV 4 , according to the evolution the time for the little rip is, 

t LR ~ 10 13 Gyrs , (69) 

which is a large period compared with the current age of the universe. For other kind of expansions, as the an 
exponential Hubble parameter (studied in pj}); this time can be much shorter (~ 300Gyrs). However, in an ekpyrotic 
scenario the occurrence of a Little Rip will depend on the duration of the accelerating phase before this ends, and 
a new contracting phase starts again. Note also that close to the dissolution of the bound structure, gravity will be 
very strong, and the breaking of Lorentz invariance will be present, such that the values of (A, n) will determine the 
expansion rate, and for instance the occurrence of the Little Rip. 

Let us now consider a model that may reproduce a entire cycle of an ekpyrotic universe, 

H = H Q - H ie - 0t . (70) 

For H\ > Hq, the Hubble parameter (|70p represents a universe that crosses through out a contracting phase, and 
then ends in an accelerating expansion for large times. Obviously, one would need to provide the way to start a cycle 
again, however for a qualitative description, we assume here that the cycle starts again after the accelerating phase 
somehow. For the solution (1701) . we have 

R = AH 2 + 6fiH = A{Hl - 2H H l e- f}t + J ff 1 2 e~ 2,3t ) - 6nl3H ie ~ pt (71) 
where we recall that A = 3(1 — 3 A + 6/u). From (fTTj) we get 



(AH + 3/x/3) ± J(AH + 3 M /3) 2 - (AH 2 - R) 
e pi = (72) 



Hi 

For simplicity we consider the case when AH + 3/x/3 = 0. Then Eq. (IT2"j) gives 



JR - AH 2 

-» = ± V = • (73) 
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And the Hubble parameter (|70|) can be rewritten in terms of the scalar curvature R, 

H = Hq — H ie -' jt = H oT y/R-AHl (74) 
In this case the first Friedmann equation (|18p yields, 

\2^(i{AH 2 - R)(H T ^R-AH 2 )f" - BF + F - k 2 Pm = , (75) 

where B = 6[(1 - 3A + 3(i)H 2 + f iH}. Then, by setting /3 = 2g ° (1 ~ 3A+3/j) , we obtain, 

B = 6(1 - 3A + 3//)(l - A)H 2 + 6(1 - 3A + 3fj,)R . (76) 

Eq. (ff5)) is still a very difficult expression, so that the search of exact solutions for F(R) is a difficult task. Nevertheless, 
we can reconstruct some particular exact actions by considering special matter fluids. Let us consider the matter 
energy density, 

Pm = k- 2 [12^{AHI - R)(H T ^R-AH 2 )F" - Ca~ 3 } . (77) 
Then the FLRW equation ([75]) admits the following particular solution 

F(R) = d [6(1 - 3A + 3fx)R + 6(1 - 3A + 3/x)(l - A)H 2 ] rn^+m . (78) 

In a similar way, other particular solutions of the Friedmann equations can be reconstructed. Hence, we have shown 
here that ekpyrotic universes can be well described in the frame of Hof ava-Lifshitz gravity. 



V. DISCUSSIONS 



In the present paper, we have analyzed some particular cosmological solutions in the context of Hofava-Lifshitz 
gravity, where basically some generalizations of the original action [8], similar to standard F(R) gravity, have been 
studied. It is well known that for a particular Hubble parameter, the corresponding action can be reconstructed in 
the framework of F(R) Hofava-Lifshitz gravity (see Ref [la]), where the presence of the set of parameters {A,/i}, 
consequence of the restriction of the symmetries of the theory, can vary along the cosmological evolution, since 
their value depends on the energy scale of a particular system. Hence, the presence of this set of parameters will 
fluctuate along the universe evolution, affecting the corresponding cosmic solution. By assuming that GR should be 
recovered when R ~ Hq « mi ~ l0 74 GeV 4 , the parameters A = p ~ 1 during the radiation/matter dominated 

epoch and the current accelerating era, while it becomes large when R cx m^, where the quantum effects should 
become important. In this sense, the effects of Hofava-Lifshitz gravity, and specifically the extra scalar mode, may 
become important when the universe reaches stages as the Little Rip, or other phases from a typical ekpyrotic universe 

Hence, in the particular solutions studied here, the ekpyrotic scenario becomes an important focus for analyzing 
Hofava-Lifshitz gravity, as the universe owns a periodic behavior, crossing different stages, where the quantum nature 
of the theory may be relevant. Moreover, we have shown that particular actions which lead to a cyclic nature of the 
Hubble parameter can be reconstructed. Several techniques have been used for the reconstruction procedure. By 
using an auxiliary scalar field, coming from the F(R) sector, we have shown that cosmological solutions can be easily 
obtained. In addition, we have studied the shape of the action along each phase of a typical ekpyrotic universe, where 
the corresponding actions have been obtained. It is straightforward to show that such actions lead to standard F(R) 
gravity when A = /i = 1 , and can be identified with some particular viable theories ■ Then, we can conclude that 
this class of actions can perfectly describe the entire universe evolution by means of an ekpyrotic model. Moreover, we 
have suggested the compatibility between an ekpyrotic universe and the presence of a Little Rip, a non singular point 
that may lead to the break of some bounded systems, where the effects of Hof ava-Lifshitz gravity turn out important, 
and A ^ 1,/x ^ 1 . Future singularities can not be compatible with a cyclic universe unless a cure for the singularity 
is considered A next step should be to probe the possibility to reproduce cyclic cosmologies within the frame of 
so-called viable F(R) gravities (see for instance, Ref. [28]). While the violation of Newtonian law can be avoided in 
F(R) Hofava-Lifshitz gravity (see [l6j])> the presence of instabilities and other features should be studied in more detail. 
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On the other hand, in order to have a complete picture of the universe evolution, one should specify how reheating 
occurs. Nevertheless, this is beyond of the scope of this work, but an interesting proposal for a reheating mechanism 
in the frame of UV complete theory is pointed out in [29( . 

Therefore, in an ekpyrotic universe, the main implications of F(R) Hofava-Lifshitz gravity would come during those 
phases when the full diffcomorphisms are broken, basically during the early and ending phases, that may affect other 
classical eras, specially by the perturbations, which should be an important point to be studied in the future, where 
the effects may be distinguishable from other models. 
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